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The ground state ensemble of the highly frustrated pyrochlore-lattice antiferromagnet can be 
mapped to a coarse-grained "polarization" field satisfying a zero-divergence condition From this it 
follows that the correlations of this field, as well as the actual spin correlations, decay with separation 
like a dipole-dipole interaction (1/|_R|^). Furthermore, a lattice version of the derivation gives an 
approximate formula for spin correlations, with several features that agree well with simulations 
and neutron-diffraction measurements of diffuse scattering, in particular the pinch-point (pseudo- 
dipolar) singularities at reciprocal lattice vectors. This system is compared to others in which 
constraints also imply diffraction singularities, and other possible applications of the coarse-grained 
polarization are discussed. 

PACS numbers: 



I. INTRODUCTION 

Highly frustrated antiferromagnets are characterized 
by a very large number of essentially degenerate states, 
such that in a range of temperatures much smaller than 
the spin interaction scale they have strong local correla- 
tions, yet fail to order 0: diffuse scattering is the obvious 
probe of such a state. In this paper, I make the point that 
the classical ground-state ensemble can entail constraints 
that ensure a generic power-law decay of correlations: a 
picture of theses state as "liquid-like" (i.e. featureless) is 
thus imcomplete. 

We consider specifically the pyrochlore lattice, con- 
sisting of corner-sharing tetrahedra, partly because of 
its simplicity as a model [high (cubic) symmetry and 
nearest- neighbor interactions] , but most importantly be- 
cause it is the magnetic lattice in many real systems: the 
B-site spinels 2, 3J, two families of pyrochlores repre- 
sented by CsNiCrFe HQ, and the oxides YzMnzOr 0, 
and even the transition-metal lattice in "frustrated itin- 
erant" antiferro- or ferrimagnetic Laves phases, such as 
Yi_xScxMn2 0,0. Ferroelectric degrees of freedom (wa- 
ter ice 0,^3) ^■'^d charge order (in magnetite 0) are 
also equivalent to pyrochlore systems. 

Different pyrochlores show a wide variety of magnetic 
behaviors at low temperatures: apart from long-range or- 
der, these include lattice distortions [2>E|) "spin-ice" be- 
havior 

El 111 El 111 113 > paramagnetism (in TbaTiaOy) 
down to low temperatures, E3 glass freezing (in 

Y2M02O7 even when structural disorder is quite small [g. 

However, most of these have a cooperative paramagnet 
regime at higher temperatures ~ but still quite low com- 
pared to the interaction scale J - demonstrating that 
they have a large set of nearly degenerate (and accessi- 
ble) states. The present paper addresses this phase, as it 
would be extrapolated to zero temperature. 

In the remainder of this section, I review the lattice and 
the ground states of a pyrochlore antiferromagnet, and 
survey three situations in which it can be modeled by 
the ensembles used in this paper. Then (Sec^jthe map- 



ping of the low-temperature states of the Ising model to 
the "diamond ice" model is to construct a coarse-grained 
polarization field, which functions somewhat like an or- 
der parameter for this model, and this is used to predict 
some unusual features of the magnetic diffuse scatter- 
ing. Sec. mil presents another version of the derivation, 
which not only gets the long-wavelength behavior (cor- 
responding to the neighborhoods of special points in the 
Brillouin zone) but provides an approximation for the 
entire zone. I also discuss (Sec. IIVII other situations in 
which local constraints in a ground state produce diffrac- 
tion singularities, and speculate on other questions that 
can be addressed using the insight that the polarization 
field contains the important long-wavelength degrees of 
freedom. 

A parallel paper contains many of the same results. 



A. Pyrochlore lattice and Hamiltonian 

The pyrochlore lattice consists of corner-sharing tetra- 
hedra arranged in cubic (fee) symmetry. Its key property 
is that the sites are the bond midpoints of a diamond lat- 
tice. [See F ig.[Tl other figures emphasizing the tetrahedra 
are in Refs.ll2lor lT^l20l ] For future use, we define vectors 
pointing towards the four corners of a tetrahedron. 
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cyclic perms. 



(1.1) 



In the diamond lattice, every even-to-odd bond vector is 
a u„i. 

The basic results derived in this paper are transfer- 
able to other lattices in which the spins sit on the bonds 
of a bipartite lattice, and with antiferromagnetic inter- 
actions among all spins on bonds sharing a common 
endpoint,such that their sum is constrained to a fixed 
value. Besides the pyrochlore lattice, this class includes 
the kagome lattice, the two-dimensional "checkerboard" 
lattice. El 113 , the "sandwich" lattice modeling SCGO 
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and consisting of two kagome layers linked by an addi- 
tional triangular layer of spins, [23, 113,01 and the garnet 
lattice j25l |. 

The spin Hamiltonian, having only nearest-neighbor 
antiferromagnetic exchange, can be cast as a sum of 
squares: 

Hspin = ^spin ^ Si ■ Sj + COnst = ^^y^ ^ L^, (1.2) 

(ij) a 

where the total spin of four sites in a tetrahedron is 

La EE ^ Si (1.3) 

where a labels each diamond site and "i € a" runs over 
the surrounding four spins. For simplicity, I will mostly 
treat the case with Ising spins ti. (There is little loss of 
generality, in view of Sec. lIBI below.) 

H = iJ^L^ = J^tjtj + const. (1.4) 

" (ij) 

From p.4|) it is evident that any state with 

i„ = 0. (1.5) 

is a classical ground state; these are massively degener- 
ate. The frustration of the pyrochlore Ising model, and 
the extensive ground state entropy implied by Ijl.SII , were 
recognized very early 

B. Temperature regimes and validity of Ising 
model 

Here I discuss the three different situations in which re- 
alistic models can be treated under this paper's scheme. 
[Readers who are content with a treatment of an Ising 
toy model, and do not immediately demand any connec- 
tion to experiment, may skip to Sec. ITU where the central 
result is derived.] 

First, it is possible that spin anisotropies (additional 
to (|1.2|) ') reduce the ground state manifold to the Ising 
states. If the lattice is to remain cubic, the preferred axis 
cannot be the same for each spin, but instead is Um(i), 
where = Aum/^/Sis the unit vector of 1)1. and m{i) 
labels the direction of the diamond-lattice bond on which 
spin i sits. Assuming the spins are strongly aligned, one 
puts Si — tiUjn(i) and the Hamiltonian reduces to (|1.4() 
with J = — Jspin/3, since u,„ • u„i' ~ —1/3 for ra ^ m! , 
and m[i) ^ 'rn{j) for nearest neighbors. This is the "spin- 
ice" model ISiMj- 

The other two situations correspond to isotropic 
Heisenberg models in the large-S* semiclassical regime 
(where S is the spin length). This is necessary to guar- 
antee that (|1.5|l is a good starting point for specifying 
the ground states. One can identify a succession of en- 
ergy scales in the pure nearest-neighbor Heisenberg anti- 
ferromagnet, corresonding to successive breakings of the 
grond state degeneracy as more effects are considered. 



The largest scale is Ej ^ zJS"^ where z is the coordi- 
nation number; this is the scale of the mean-field (class- 
cial) part of the exchange energy, and of the Curie- Weiss 
constant. At temperatures T ^ Ej, it is reasonable to 
approximate the spin ensemble by a subset of the ground- 
state ensemble. 

The spin- wave spectrum {ojj} depends upon which 
state we linearized around. Thus, the total spin-wave 
zero-point energy (to harmonic order, "Y^- \h,ujj) partially 
breaks the gro und-state degeneracy, favoring coUinear 
states, 0,f2g(. This is expressed quantitatively by an 
effective Hamiltonian j22L.27^.3Q.] . with coupling constant 
~ E'coii, where the next largest scale is EcoW ^ JS Ej. 
(Even when S is not so large, e.g. 5" = 5/2 as with real 
spins, the coefficient of JS is commonly less than 1/10 so 
the inequality is valid.) In lattices other than pyrochlore, 
one finds a similar regime in which a discrete subset of 
the classical ground states (e.g. coplanar states, |^ 0| 
in the kagome case) gets selected by the harmonic zero- 
point energy. 

Finer treatment of the spin- wave fluctuations produces 
an even smaller energy scale iodise for selection among the 
discrete states. In the kagome case, [s^l^; this is due to 
anharmonic terms in the Holstein-Primakoff expansion, 
so that iodise c>c S^^^ <^ S. In the pyrochlore case, besides 
the anharmonic terms, E^isc does include contributions 
of 0{JS) from the harmonic terms. ■ However, these 
terms only partly break the degeneracy, and the coeffi- 
cients are small, so we can still assume -Edisc ^ -E'coU- 

Thus, the second regime in which the Ising theory mod- 
els the correlations is iodise < T < i?coii- One expects 
a symmetry breaking to long-range coUinear order, in 
which a global axis n is spontaneously adopted and each 
spin is given by = t^n, with ti = ±1, plus small fluc- 
tuations. This phase is a kind of "spin nematic" [s^ . 
The energy terms distinguishing different Ising states are 
unimportant since T > i?disc, hence the spin ensemble is 
roughly the ground states of H1.4|l . as claimed. 

Of course, in the real world there is another en- 
ergy scale -Epcrt for perturbations of the pure Heisen- 
berg Hamiltonian, representing the magnitude further- 
neighbor exchange and dipole couplings, as well as lat- 
tice distortions, and disorder, which cause a transition 
or freezing into some other state at sufficiently low tem- 
perature, so the regime of the spin nematic-like phase is 
really max(£'disc, Export) <T < E^oU- 

Finally, the third situation is when -BcoU ^ T <C Ej. 
This ensemble is well modeled by the low-T limit of the 
classical Heisenberg model. It turns out (Sec. IIIDI be- 
low) that the key notion of "polarization" does extend to 
the Heisenberg spin case. I have mostly neglected that 
case in this paper because, at my level of approximation, 
the results look identical; it would merely obscure the 
notation by adding spin-component indices everywhere. 
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II. COARSE-GRAINING, FOURIER MODE 
FLUCTUATIONS, AND LONG-RANGE 
CORRELATIONS 

In this section, I present the steps leading to power-law 
correlations using the framework of a continuum theory, 
in the ideas appear more transparently. 



A. Ice mapping and local polarization 

In fact the pyrochlore (Ising) ground states map 1- 
to-1 onto those of the well-known diamond-lattice ice 
model 0, [S^ , in which the degrees of freedom are ar- 
rows along the lattice bonds. In the map, every tetra- 
hedron center becomes a vertex of the diamond lattice, 
while the spin sites map to bond centers of the dia- 
mond lattice. Each spin ti = -|-1(— 1) maps to an ar- 
row pointing along the corresponding diamond lattice 
edge, in the positive (negative) sense from the even to 
the odd vertex. This well-known mapping pj) . Issl Is^ 
is also used to model the "spin ice" system Dy2Ti207 
(and also Ho2Ti207), wherein local (111) anisotropy 
plus ferromaqnetism makes a highly frustrated Ising 
model ,13„ 1^ J^J-) The ground state condition - 
net spin of every tetrahedron is zero - maps to the "ice 
rule" the numbers of incoming and outgoing arrows are 
equal at every vertex |^. 

The key object in this paper is the ice polarization field. 
On a diamond-lattice bond a polarization tiU„j(i) can be 
defined, aligned from the even to odd diamond-lattice 
vertex if the spin is up, oppositely if it is down; here 
m(i) is the local 3-fold axis of site i. On every diamond 
vertex, we define the mean of this polarization over the 
surrounding tetrahedron of spins: 



P(R„) = ^t,u„ 



(2.1) 



The six possible ground-state configurations of that 
tetrahedron correspond to P(Rq) — (±1, 0, 0), (0, ±1, 0), 
or (0,0, ±1). 

Finally, the coarse-grained arrow field P(r) is the 
coarse-grained polarization averaged over some larger 
neighborhood and assumed to vary smoothly. 



B. Effective free energy and correlations 

The ground-state entropy density is a function of the 
average polarization. A subensemble of states in which P 
is large (which can be forced by boundary conditions) has 
relatively little freedom for rearrangements of the spins 
or arrows; indeed, for a saturated polarization such as 
P — (1, 0, 0) the ensemble consists of a single microstate. 
Thus it is very plausible that the entropy density has a 
maximum for zero polarization. Therefore, to lowest or- 
der, the total free energy (arising entirely from entropy), 



as a function of coarse-grained P(r), has the form 
Ftot({P(r)})/T = v^^lj d'r^K\P{r)\\ 



(2.2) 



where Uceii = a"^/4 is the volume of a primitive unit cell. 
The "stiffness" k is dimensionless, as appropriate since it 
is purely entropic in origin. 

Corresponding to the condition (|1.5|l . i.e. the ice rule, 
P(r) satisfies a divergence constraint 



V • P(r) = 



(2.3) 



like a magnetic field without monopoles. Eqs. H2.2|l and 
(|2.3|) look, respectively, like the field energy of a mag- 
netic (or electric) field, and its divergence constraint, in 
the absence of monopoles (or charges). These equations, 
together, signify that that the probability distribution of 
the (long-wavelength portion of) the polarization field is 
the (constrained) Gaussian distribution 



Prob({P(r)}) (X e 



-F({P(r)})/T 



n'5(V.P(r)) (2.4) 



Fourier transforming H2.2|l simply gives Ftot = 
X^k a naive use of equipartition would give 

(P^(— k)P,y(k)) = (1/k)(5^i,. But the divergence con- 
straint (|2.3|) imposes 



k-P(k) = 



(2.5) 



in Fourier space. Thus the correct result has the longi- 
tudinal fluctuations projected out: 



(P^(-k)P.(k)) = i(. 



(2.6) 



Fourier transforming H2.6(l back to direct space gives 

(2.7) 



47r r 1 
(P^(0)P,(r)) ^ — 6{r) + —{Sf,, - 3f^f,) 



at large separations r (where r = r/|r|.) Correlations 
have the spatial dependence of a dipole-dipole interac- 
tion, which is a power law. Models that exhibit such 
correlations - including the pyrochlore system, so often 
described as "liquid-like" , - are thus, in a sense, in a crit- 
ical state. [The generalization of (|2.7|) for d-dimensional 
real space, would be a l/r"^ decay.] 

This criticality was first appreciated in the ice model 
itself, being detected originally in a simulation [s^- A 
functional form with a dipolar singularity like II2.6I) was 
produced by a clever random-walk approximation to a 
series expansion |39j . The universal explanation, made 
here, that dipolar correlations arise from H2.2(l with the 
divergence condition, was first put forward to explain 
experiments on two-dimensional ice- like systems |l0(. 

Ref. have also presented an ansatz equivalent to 
(|2.4|) . derived the dipolar correlations, and confirmed 
them by simulations, for the dimer covering of a sim- 
ple cubic lattice. [See their eq. (3).] The l/|r|'^ decay 
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has also been obtained analytically and numerically in 
Ref. l6Q | ■ f or the pyrochlore model of this present paper. 

Ref. 12(1 Sec. II D 1, recognized that the ground-state 
constraint in the Heisenberg pyrochlore antiferromagnet 
(|1.5() entails long-range correlations, but in the absence of 
the polarization concept, the argument's form is diffuse, 
and it was not possible to predict an explicit functional 
form. An interesting heuristic argument was made there 
to justify the empirical fact of the "bowtie" shape (see 
Fig. |2Jl taken by the diffuse scattering: in other words, 
that it takes a scaling form in terms of g^/gj.. 

C. Diffraction consequences 

1. Spin structure factor 

There is a linear relationship between the actual Ising 
spins and the ice arrows: specifically, the map H2.1|l . from 
spins to the the diamond-vertex polarization is actually 
invertible: 

t(R±iu„,)-4u,„-P(R), (2.8) 

where we take the -I- or — signs for even and odd ver- 
tices R, respectively. Thus it is not surprising that the 
pseudodipolar correlations H2.7|l of the latter imply sim- 
ilar 1/r^ correlations for the former. However, the co- 
efficients relating these (namely, the {u^} vectors) are 
staggered in sign. 

To see this, first Fourier transform (|2.8|) . using the fact 
that 

u„^. = ie5'^--("™-"-\ (2.9) 

similarly uu^^.y with K020, etc. When this is inserted 
into the formula for i(k) [similar to 1)3. l|l ]. the staggering 
contained in the u,„ factors shifts the argument by the 
reciprocal space vector K200 so that i(K2oo+q) -fx(q), 
etc. Consequently, the singularities of the spin structure 
factor in reciprocal space are pseudo-dipolar in form, just 
like (|2.10() , but occur at nonzero reciprocal lattice vectors 
K rather than at q = 0. 

Let i(k) be the Fourier transform of the spins. The 
structure factor (i.e. the neutron diffraction intensity, 
modulo polarization factors) has pseudo-dipolar singu- 
larities 

(|i(K2oo + q)P) « (2.10) 

^ll ' 9_L 

Here K200 = [27r/a](2, 0, 0) is a reciprocal lattice vector; 
and q|| = and q\ = qy + q'^. Near Km = 27r(l, 1, 1), 
the same form holds with qy and qj_ being the compo- 
nents of q parallel and perpendicular to the (111) direc- 
tion. 

The elastic constant k must be determined from Monte 
Carlo simulations. Yet, without knowing it, we can still 
make the nontrivial prediction that the diffuse scattering 



has quantitative the same strength near Km as it does 
near K200 • In an isotropic Heisenberg model, these pseu- 
dodipolar singularities can be distinguished from true 
dipolar singularities in neutron scattering experiments, 
since they are independent of spin direction. They can 
be seen in only one spatial direction around each recip- 
rocal lattice point. 

The functional form (|2.1U|) has nodes which we expect 
(correctly) are a consequence of symmetry and therefore 
extend throughout reciprocal space, beyond the small-q 
limit in which H2.1(J|I was derived. For example, ()2.1(J|I 
indicates that 

(|f(A:„0,0)|2) =0 (2.11) 

for wavevectors anywhere along the entire (100) axis. 
Similarly the structure factor is zero along the (111) axis 
near Km from (|2.10|l Since k = (000) is the intersection 
of seven distinct axes, along each of which the structure 
factor must vanish, we deduce that the diffuse scattering 
vanishes near the origin as 

(|t(k)|2) cx {klk^ + 5 permutations) - Qklklkl. (2.12) 

Indeed, experiments and all simulations on pyrochlore 
antiferromagnets find quite small diffraction throughout 
the first Brillouin zone. 



2. Derivation of nodal lines 

It is possible to confirm H2.11|l without using coarse- 
graining, by sharpening an argument made in Sec. II D 
of Ref. .201 (I have written the steps in detail, so it will be 
clear whether they do or do not carry through for other 
lattices.) 

Partition the pyrochlore lattice sites into (100) planes 
labeled by x, and define tpiano(2;) = J2x =x^i- Then 
consider, for example, the tetrahedra with centers at 
a; = 0: each such tetrahedron includes a pair of spins 
with Xi = — J and another pair with Xi = +j. By the 
ground-state constraint 11.5|l the sum of one pair is the 
negative of the sum of the other pair. Since every spin 
with = ±i belongs to exactly one such tetrahedron, 
this implies tpiane(-|) = -ipiano(|); and by Bravais lat- 
tice periodicity, ipianc(i + ^n) = (-l)^<pianc(|)- Since 
all spins have Xi = j + for some n, 

i(fc,,0,0) cx ^e-''=-(i+5")(-l)"tpi,,,(i) (2.13) 

n 

which cancels, except for a possible Bragg peak when kx 
is a multiple of 2tt. However, in the maximum entropy 
state that we assumed in Sec. IIIBI the average spin in 
one of these planes - which is proportional to the mean 
polarization - is zero, so i{kx, 0,0) = for all k^- 

A similar, but more complicated, argument works for a 
{111} direction. Let = • :^(1, 1, 1) be the projection 
of site i on the (111) axis, and redefine ipianc(C) as the 
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spin sum over the spin plane with = ^. Each tetrahe- 
dron with ^ = includes three spins in the ^ = — 1/ 4\/3 
plane and one spin in the ^ = \/3/4 plane, and every spin 
in either plane belongs to a unique tetrahedron; conse- 
quently tplanc(-l/4V3) = -tplane( + %/3/4), CVCU thoUgh 

one plane has three times as many spins as the other. 
Furthermore these unequal planes are equally spaced, so 
ipiane(\/3/4 + n/VS) = (-l)"ipiane(y3/4) and the rest 
of the argument follows as before. 

3. Alternate derivation of nodal lines 



replacement 7i — > — H • is equivalent to substi- 

tuting 

Lq, — > Lq, — H/2Jspin (2-17) 

in (|1.4|l . In the case of vector spins (but not Ising spins!), 
one achieves a ground state by satisfying Lq, = H/2Jspin 
on every tetrahedron. Simply replacing 6si = — H/8 
in the definition of P(r) (i.e subtracting off the mean 
spin expectation), provides a polarization appropriate to 
this ground state, which should exhibit the same sort of 
long-range correlations (but with an |H| dependence.) 



There is an alternate way to see that diffuse intensity 
must be small near the zone center. Let P"*" (k) be Fourier 
transform of the discrete polarization on (say) even di- 
amond sites, defined like H3.2|l . On the one hand, the 
Fourier transform of the constraint H1.5|l at odd vertices 

4G(-k)-P(k)^0 (2.14) 

where 

4 

G(k) = ^ e*'^ "™u„ « i*k (2.15) 

m— 1 

at small wavevectors. On the other hand, we can start 
from (|2.8|) and take its Fourier transform - this time, not 
taking advantage of H2.9|l : the result is 

i(k) =4G(ik)-P(k). (2.16) 

Near k = 0, comparing H2.16|l and (|2.14|) and recalling 
that G(— ik) w — iG(k), we conclude that i(k) w 0, as 
I asserted earlier. Furthermore, whenever k is along a 
{100} or a {111} symmetry axis, G(k) is parallel to that 
same axis by symmetry (for all k) , hence structure factor 
(|2.11l) is null all along those axes. 

D. Generalization to Heisenberg spins 

In the case of isotropic n — S component spins, we 
can reiterate all the coarse-graining arguments of Sec. ^ 
in terms of the Heisenberg spin components. Polariza- 
tion components are defined as in (|2.1|) but now for each 
Cartesian spin component, so the polarization field is now 
a tensor carrying not only space but also spin indices: 
(Pa)/3M = Ejgq(^«)/9("™(^))m- As in the Ising case, it is 
easy to convince oneself that the entropy density is max- 
imum when P = 0. Hence we expect that (|2.2|l remains 
valid, except that each |P(r)p is now interpreted as a 
tensor norm (sum of squares of all tensor elements). Fi- 
nally, (|2.3|l now becomes three equations, one for each 
spin-space component. 

The basis of this paper - and the signatures implied 
in the diffuse diffraction - extends even to Heisenberg 
pyrochlore antiferromagnets in a magnetic field H. The 



III. LATTICE APPROXIMATION 

The continuum theory presented in Sec. ^] can only 
predict the shape of singularities at special points in re- 
ciprocal space. For a better comparison to experiments, 
a theory of the diffuse scattering throughout reciprocal 
space is desirable. Of course, whereas the result (|2.7f) 
is universal across a class of models as enumerated in 
Sec. Ill Al and Sec. IIV (TTl the detailed shape of the scat- 
tering calculated here is specific to the pyrochlore lattice. 

The derivation depends on lattice Fourier transforms. 
This unfortunately forces the introduction, for this sec- 
tion, a new indexing t(R-|- ^u^) and i(R) for the same 
objects ti and defined previously, where R designate 
the diamond sites; the even sites R form an fee Bravais 
lattice and we let N be the number of primitive unit 
cells. Thus each pyrochlore site is written R -I- ^u^; 
the surrounding odd diamond vertices are at R -f u,„ . 

i^(k) = 4=J2 e-'''-^^+^''"'^t{R + iu„). (3.1) 

^ R 

The definition (O is rewritten L(R) = E1=i ± 
^u„j)) taking -I- and — when R is an even or odd diamond 
vertex, respectively. Thus the Fourier transform of (|1.3|l , 
restricted to even or odd vertices as indicated by the "±" , 
is 

L±(k) = ^^e±^'''""f™ (3.2) 

V m=l 

A. Derivation of fluctuations 

The simplest way to describe the approximation is first 
to imagine a distribution of spins 

Prob({t,}) oc e ^'o exp(-/3effH), (3.3) 

where H is given by H1.4|l . and f3cS ought to be the in- 
verse temperature. The Gaussian factor in (|3.4|l may be 
viewed, in the Bayesian spirit of the maximum- likelihood 
approach, as a trivial a priori independent distribution 
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before we account for any spin interaction. The sec- 
ond factor enforces spin correlations; here is the 
Hamiltonian H1.4|l for Ising spins, except that now ti is 
aUowcd to take any real value. 

Adopting the limit /3off oo, H3.3|l reduces to 



Prob({tO) 



■T ^ 



(3.4) 



The second factor in (|3.4() imposes the ground state con- 
straints (|1.5|) around every tetrahedron, both even and 
odd. 



1. Ising spin correlations 

Our goal is to evaluate {Utj), and this is not hard 
to do using the Fourier transform (13.11) . We must first 
rewrite Uc^SiL^) cx ll^^S{L+{k))S{L- (k)) in and 
re-express this in terms of {tm(k)}. 

Since different wavevectors decouple, it will be conve- 
nient to view tmO^) as a complex 4-component vector. 
Then the even and odd vertex constraints take the form 
of orthogonality conditions: 



L±(k) = J2 e±2''""'"i;„(k) = (£;±(k,t(k)) = (3.5) 

m— 1 

The coefficients are 

£±(k) ^ e^r'--™ (3.6) 
so = (£^^)*. The inner product is defined by 



B„ 



(3.7) 



m— 1 



So, the distribution H3.4|l can be rewritten in Fourier 
space as 



n 



S{{E+(k),i)S{{E-(k),i) 



(3.8) 



Thus, for each k in (|3.8|) . we now have a Gaussian dis- 
tribution over a four-dimensional space, with two con- 
straints reducing it to a two-dimensional subspace. 

Our object is to obtain correlations of the form 
(t;(— k)f,„(k)). The result will essentialUy be the pro- 
jector n into the subspace satisfying l|3.5|l . If we define 
the 4x2 matrix E, the columns of which are E'^ and 
E^ , then 4x4 projection matrix for the (desired) space 
orthogonal to E+ and is H = / - E{E^E)-^E^. In 
fact 



{E^E) = 4 



{E^E)-^ = -{l-\H\ 



1 

H* 
1 

-H 



H 
1 

-H* 
1 



(3.9) 
(3.10) 



where the k argument in E and H was suppressed; here 

= cos COS — COS 

4 4 4 

-i sin — sin — sin — . (3.11) 

4 4 4 



ff(k)^i^e^'^- 



[H{-k) = H{k)*, |i?(k)|2 ^ 1(1 + cos^cos% + 
cos ^ cos ^ + COS ^ COS ^). The final result is 

(t,(-k)t„,(k)) = 4 {Si„, - [EiE^E)-^E%n} (3.12) 

where E and {E^ E)^^ are implicitly functions of k and 
are defined by (EH) and ^J^. 

Since Trll = 2, wc must have = 2 to satisfy the 
normalization condition (tf) = 1. 



2. Spin structure factor 

The structure factor, as measured in neutron diffrac- 
tion, combines the four sublattice contributions. Write 
the Fourier transform of all the spins as t{k) = 

^E™ERe-^'^-^''+^""'^(R + K") = {M,i{k)), 

where we used the 4-vector inner product and 



M = (1,1,1,1) 



(3.13) 



The structure factor is then a projection of (|3.12|) onto 
the M vector, thus 



(|i(k)|2) = tlM^[I - E{E^E)-^E^]M. 



(3.14) 



In the case of spin ice, the mapping of the Ising spins 
{ti\ to the real spins is different (it is staggered). In 
this case, Mpm — {'^m)f3- The pinch-point singularities 
will appear, but at different places (including the origin), 
because in this case the spins are the local polarizations. 
In this case, the scattering around the origin is no longer 
suppressed. 

In either the antiferromagnet or the spin-ice case, there 
will be extinctions of the expected singularity at points 
where {E^,M) happens to cancel. 

Considering the form of H3.12|l , singularities in the spin 
fluctuations may (but do not necessarily) occur when the 
2x2 matrix (E^E) is singular. In view of (|3.10() . this oc- 
curs when |i/(k)p — 1. But that is precisely the defining 
condition for the fee reciprocal lattice vectors K. [This 
was derived by a different route in Sec. Ill C II ] 

Indeed, £;+(K2oo) = i(-l, -1,1,1) and i;+(Kiii) = 
^(—1, 1, 1, 1). We see E~^(K) is proportional to a real 

four- vector. [That generically happens at a point in 
three-dimensional k-space, since there are three indepen- 
dent phase relationships to be satisfied among the com- 
ponents {£'+}.] Then i?+ is proportional to E^ , and the 
rank of i?(K) is reduced from two to one at these points, 
confirming that E(K) is singular. 



7 



The special hnes in reciprocal space on which diffrac- 
tion is zero, are those where the 4- vector of coefficients 
(13.13(1 . relating the physical spin to the four sublattice 
spins t„i(k), happens to be orthogonal to both of the 
vectors in the null space of the 4x4 matrix H3.12|l . 

B. Finite temperature 

In the approximations of Sec. HTIand Sec. IIII^ where 
the tetrahedron constraint was imposed rigorously, we 
were forced to project out the corresponding fluctuation 
modes. It is possible to extend the approximation so as 
to permit small fluctuations of the tetrahedron magne- 
tizations, as must be excited at T > 0, starting from 
(|3.3() . Now, no projection is required, since a large co- 
efficient PcS tends to suppress those fluctuations, and 
so this PcS < oo case is actuallly more straightforward. 
[This subsection is essentially a streamlining of Ref . Isol 
as converted to my notations.] 

The Hamiltonian in (|3.3|) is a quadratic form in the 
ii's. Fourier transforming (| 1 . 4|l . and putting it into the 
4- vector notation with (|3.5|) . we get 

n = yJ2 + |L-(k)p] = lj{i,E^Ei). 

(3-15) 

Thus, 1(3.3(1 becomes Prob({im}) oc exp(— ^(i. At), where 

A=^I + f3,sE^E. (3.16) 

and 1(3.12(1 gets replaced by 

(iz(-k)i„(k)) = {A-^)i.m. (3.17) 

The matrix E^E has rank two. It can be seen that, as 
PeS — > OO, ((3.17(1 indeed reduces to the projector ((3.121) . 

The consequence of /3off < oo for correlations is that, 
in reciprocal space the pseudodipolar singularities ((2.10(1 , 
get rounded by the substitution 

q^^q^ + r^ (3.18) 

1 /9 

in the denominator, where ^ ~ PcS is a a correla- 
tion length. In real space, the power-law decays ((2.7(1 
acquire an exp(— r/^) factor that cuts them off, as noted 
in Ref.El 

C. Other analytic approximations 

Three prior treatments of the diffuse scattering in the 
pyrochlore arrived at a mathematical form more or less 
identical to ((3.3(1 . but with different formulas for the co- 
efficients in these equations as a function if temperature. 
Thus, of course, their result is ((3.17(1 : however, they did 
not note the pseudodipolar (or, in the spin-ice case, lit- 
erally dipolar) correlations which are implicit in these 
formulas at the T = limit. 



The diffuse scattering problem was first addressed by 
Reimers [4^ for Heisenberg, or in general n-component 
vector spins.) That derivation is based on mean field 
theory (Ornstein-Zernike correlations), which ought to 
be valid in the critical regime T , where — 

J/n. However, that is invalid in the pyrochlore case, since 
the real critical temperature is driven to zero. Even so, 
the k dependence of the result (eq. (15) of Ref. ,49), for 
the hmit T T^^ , is exactly of the form with 
1/4 = 3(r/Tf F - 1) and /Jeff = 1/T. 

Canals and Garanin |47l li^ considered the classical 
pyrochlore antiferromagnet with n-component unit spins, 
in the large-n limit which is tractable analytically. This 
is, in effect, like the "spherical model" approximation for 
finite-n, in that a constraint |si| = 1 on every spin is 
replaced by one on all N spins, ^ . js^p = N. As they 
noted in Ref. 0, Sec. Ill, the n = oo limit is completely 
described by Gaussian approximations. Their result for 
small temperatures can be reduced to (I3.3() with I/Iq — 
n/3 and PcS = 1/T; this reduces to (|^ as T ^ 0. 

Finally, Yoshida et al [s^l use an elaborate cluster- 
variational method to derive a sensible formula for the 
temperature-dependent diffuse scattering for the Ising 
in = 1) pyrochlore antiferromagnet (specifically, spin 
ice). Their result, taking the lowest-order finite tem- 
perature correction, amounts to ((3.3(1 with = 2 and 
l^csJ = i ^e^'^ ^'^, which in the zero-temperature limit re- 
duces to [Note that Joff/3 of Ref. HO is my J. My 
i?*™ are essentially linear combinations of the 4-vectors 
c and s of their eq. (B-3).] 

Until now, I have only mentioned approximations 
which involve the connectivity of the sites and which 
can incorporate the long-range nature of the constraint. 
Some other approximations, intended mainly to model 
the bulk susceptibility are based on a single isolated 
tetrahedron [S^, |52, |53 ■ Apart from the restriction that 
this tetrahedron has zero net spin, these approaches nec- 
essarily miss the long-range constraint and hence give a 
poor picture of the long-range correlations. 



D. Comparison to diffraction in experiment and 
simulation 

Several pyrochlore systems show the same characteris- 
tic diffuse scattering features: (i) the entire first Brillouin 
zone has a very low intensity (ii) intensity is zero along 
{100} abd {111} axes; (iii) there is a pinch-point singu- 
larity of form ((2.10() at K200 and also Km reciprocal 
lattice vectors. It should be noted that, since there are 
four spin sites per unit cell, the diffuse scattering is not 
periodic with the Brillouin zone. 

ExperimentallyjSuch features were seen in an itinerant 
Laves phase (Ref.H, Fig. 3); in the pyrochlore CsNiCrFg" 
(Ref. 1^ Fig. 4); and most recently in the spinel ZnCr04 
(Ref. I3, Fig. 3(a,b)), in the higher-temperature regime 
above a structural transition. In simulations, such pat- 
terns appeared in Fig. 2 of Ref. 0, Fig. 4 of Ref.|^ and 
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in Ref. li^. For comparison, spin-1/2 results from exact 
diaagonalization are shown in Ref. fs^. Fig. 4; they are 
quahtatively similar, but less sharp. 

Images of analytic large-n calculations of Garanin and 
Canals can also be compared: Ref. |48j. Fig. 4 (which is 
the kagome system) and Ref. |5^ Fig. 6. 

Ref.im in their Fig. 5, plot diffuse scattering from sim- 
ulation of a two-dimensional spin problem equivalent to 
the honeycomb dimer covering. (This is a plane of spins 
in the "kagome -ice" phase, whereby an external field 
applied to a "spin-ice" pyrochlore system causes {111} 
planes to decouple). The pinch-points (called "bowties" 
by those authors) are prominently visible, which are di- 
agnostic of pseudodipolar correlations on real space. 

The structure factor (|2.10() has a local maximum at 
q + X = 0, if we vary along a line q_L = const, off- 
set slightly from the g^-axis. The same behavior is found 
around Km and, of course, all other symmery- equivalent 
reciprocal lattice vectors: the structure factor has max- 
ima in the plane perpendicular to the radial direction in 
reciprocal space. The union of these planar facets forms 
the same shape (a cuboctahedron) as the boundary of 
the FCC lattice's first Brillouin boundary, but doubled 
in all three directions. Indeed, in a plane of reciprocal 
space from a single crystal of Yi_j;Sc2;Mn2, the diffuse 
scattering is concentrated near the lines where this plane 
cuts the Brillouin zone faces 8]. 

Powder diffraction data from pyrochlores showed a 
characteristic maximum pti | which was seen experi- 
mentally ij], and in simulation jl^. This is consistent 
with the fact that the (doubled) Brillouin zone bound- 
ary (where diffraction is maximum along any ray in re- 
ciprocal space) is roughly a sphere. (Note that powder 
averaging, in the vicinity of K200J amounts to integrat- 
ing (|2.1U|I over qy and q^, which yields a weakly cusped 
function Const — 27rg^ Inqx-) 

IV. DISCUSSION 

In summary, it was found that a polarization can be de- 
fined in pyrochlore antiferromagnets which - in a ground 
state - exactly satisfies a divergence condition. The 
coarse-grained version of this is analogous to an order pa- 
rameter, in being the natural variable to describe large- 
scale properties. The correlations were found (Sec. 
to have a pseudo-dipolar form which - since this is a 
pure power law - implies an infinite correlation length. 
These behaviors were repeated in a lattice-based deriva- 
tion fSec. mil In the rest of this section, I discuss other 
problems to which these findings or approaches could be 
extended or related. 



A. Implications for diffraction experiments 

The argument of this paper suggests that the analytis 
of diffraction experiments ought to focus more on the 



characteristic features, such as pinch points, identified in 
Sec. Ill CTI and Sec. IIII Deviations from the predic- 
tions at those places are sensitive measures of the extent 
to which the tetrahedron constraint is violated in the 
actual ensemble. Thus, it is suggested to analyze these 
experiments so as to extract the correlation length ^, and 
to check how well the diffraction is suppressed along the 
predicted nodal lines. [A correlation length was extracted 
in Ref. 44 from simulations, however in this case it was 
actually the finite size cutoff.] 

Deviations in the overall pattern from the shape pre- 
dicted in Sec. IIII A 21 (which are consistent with simu- 
lations '21II) are likely to reflect additional terms in the 
Hamiltonian, as in Ref.H 

Ref. j| have fitted the diffuse intensity as the Fourier 
transform, not of a single tetrahedron, but a single loop 
of 6 spins. However, contrary to the speculation in that 
paper, this does not necessarily indicate a physical state 
built from such hexagons. To explain this, I will out- 
line an alternative, equally systematic way of fitting the 
diffuse diffraction data. 

The constraint H1.4|l implies a similar constraint on the 
matrix of correlation functions. Then one can express 
any valid correlation function using a basis of linearly in- 
dependent, orthogonal functions in real space satisfying 
this constraint, in the spirit of e.g. spherical harmonics. 
The first of these terms is the same correlation that de- 
rives from the ring of six. The form observed Q is, from 
this viewpoint, the simplest possible shape, as one might 
expect at higher temperatures when all the other terms 
are damped out. To produce a pinch point, an infinite 
number of such functions would be needed, correspond- 
ing to a large spatial extent. 



B. Dynamics 

It is well known that violations of the tetrahedron con- 
straint (due to disorder, or thermal excitation) map to 
electric charges (in the language where P(r) is an "elec- 
tric field." A defect in which three arrows point outwards 
has a "charge" of Q = +1, or Q = — 1 if three point in- 
wards [S^l; in general, Q = La/2. These "charges" can 
only be created in opposite pairs; such a pair feels an 
(entropic) effective interaction which behaves (at coarse- 
grained distances) exactly like a Coulomb interaction. At 
T > 0, the interaction would be screened in the usual 
fashion, and the screening length is the correlation length 
mentioned in Subsec. IIV Al 

Apart from these defects, the total (classical) spin is 
constrained to be exactly zero. [The total magnetization 
is i J2a = by (|1.5|l .] Hence, the defects are central 
to any theory of the magnetic relaxation, as observed 
by inelastic neutron diffraction, or in AC susceptibility 
measurements [stIIss) . 

To sharpen this point, note that within the ground 
states, there is no local move that produces another valid 
state: an entire loop must be updated at once 0|. But 
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the movement of a charge along a loop leaves behind the 
same change and thus implements this nonlocal "flip" 
operation. 

The interpretation of relaxation experiments ought to 
be cast in terms of the diffusion and recombination of 
pseudo "electric charges" . The relaxation rate of the real 
magnetization is proportional to the drift mobility of the 
"charges" . The theory of their behavior is isomorphic to 
an intrinsic semiconductor; the cost a Q = ±1 defect (in 
an Ising model) is 2 J, so 4 J plays the role of the bandgap. 

Nonmagnetic impurity sites 0| act like impurity levels 
in a semiconductor (exactly at midgap). In the Ising 
model, a tetrahedron with one missing site has a ground 
state with = ±1, corresponding to Q = ±1/2: it is 
as if a quenched charge of —1/2 has been placed there, 
with the possibility of binding a +1 charge to it. In a 
Heisenberg model, such a tetrahedron often still satisfies 
Lq = and the behavior is more subtle [23| . 

We could also make predictions for the dynamic neu- 
tron structure factor. The polarization P(r) is a con- 
served quantity, so (in a classical model) it must diffuse. 
Thus 

^ = -v-Jp = rv2p (4.1) 

dt ^ ' 

Hence, near a pinch point K, this [classical] dynamics 
implies vanishing widths in the dynamic structure factor 
iS(K -(- q, w) ~ Fq^. Dynamical conclusions were derived 
by Moessner and Chalker directly from an equation of 
motion in terms of Lq, . [The consequence in a quantum 
model appears to be a gapless excitation with linear dis- 
persion |6(]|.] 



C. Other models 

1. Other three-dimensional lattices 

The polarization construction can be generalized to 
any model in which (i) "spins" sit on the edges of a bi- 
partite graph, and (ii) around every vertex, the sums 
of the "spins" is the same; here "spin" degree might be 
structural as well as magnetic. 

In three dimensions, the best examples (apart from 
the pyrochlore lattice) are antiferromagnets on the garnet 
lattice or dimer coverings on the simple cubic, bcc, 
and diamond lattices. Dimer coverings of the diamond 
lattice might realized by the ice model in which ions of 
difference valences sit on the even and odd diamond sites; 
they also correspond to the ground states of an Ising 
pyrochlore in an external field. This last system is most 
plausibly realized by a 1/4-occupancy lattice gas, with 
nearest-neighbor repulsion, representing a charge order 
problem. Alongside the dimer models are vertex models 
on the simple cubic or triangular lattice, in which each 
site has three inwards and three outwards arrows (6Qi . i61i |. 



2. Two-dimensional models 

In two dimensions, constraints such as 1)1. 5|l are en- 
countered in several models, notably in two-dimensional 
ice (= six- vertex model), the triangular Ising antifer- 
romagnet ground state, the square lattice dimer cover- 
ing, and especially the kagome Heisenberg antiferromag- 
net |6^ ISa, E3- To compare the last of these to the 
pyrochlore Ising ground states, it is fairest to consider 
the ground states of the three-state Potts model anti- 
ferromagnet on the kagome lattice. In that system, as 
in the pyrochlore, one can predict the structure factor 
(|t(k)|^) — along the lines through the origin and its 
first star of reciprocal lattice vectors; and here too the 
scattering tends to concentrate along a surface which is 
the Brillouin zone boundary, scaled up by a factor of 2. 

Following the analogy to the pyrochlore, one would 
guess the scattering has a smooth maximum at the zone 
corner Q = {|, |} type points, but in this respect point 
the behavior differs from d = 3. To understand that, 
consider that 1)2. 3|) is satisfied formally by writing 

P(r) = V X h(r) (4.2) 

where h(r) is a "vector potential". In d = 2, there is no 
gauge freedom: for a given configuration {P(r)}, h(r) is 
uniquely determined (apart from a constant), and can be 
visualized as parametrizing a (rough) interface in a 2 -I- 1 
dimensional space j^. Following standard "Kosterlitz- 
Thouless" (also known as "Coulomb-gas") prescriptions, 
spin operators have a component which is a periodic func- 
tion of the local h(r). This implies correlations with a 
power-law proportional to 1/k, and the structure fac- 
tor must have a power-law cusp at Q (which I will call a 
"zone corner sin gula rity" ) reflecting this quasi-long-range 
order [g^. Ref. l4ll have noticed these "zone-corner sin- 
gularities" in simulations; they also appear when one di- 
rectly measures the fluctuations (|/i(q)|^) of a discretely 
defined height field ^ ^ . 

In such "height models" it is also possible that the free 
energy favors a state with bounded fluctuations of h(r), 
corresponding to long-range order of the spins. 42'^ 6a|- 
In the rare cases of a height model in d = 3 (e.g. three- 
state Potts antiferromagnet on the simple cubic lattice), 
or in the ground state oi a d — 2 quantum system, this 
"locking" behavior is always expected, except in some 
quantum models which contain nontrivial Berry phases. 

Ref. have calculated the asymptotic correlation 
function for the dimer coverings of a triangular lattice, 
a system which is realized in the "spin-ice" class of py- 
rochlore in a magnetic field oriented along (111) : see their 
eqs. (4.3)-(4.4). It should be noted that two-dimensional 
dimer models (which are solvable by the methods used 
for the Ising model, i.e. free fermions) , have the peculiar- 
ity that the correlations arising from the "height" field 
h(r) have exactly the same decay (l/r^) as the pseu- 
dodipolar terms; the contributions can be distinguished 
because the first kind of correlation does not depend on 



10 



the orientation of the vector r between sites, while the 
pseudo-dipolar kind does. 

3. Quantum models 

For the 5 = 1/2 case, there is believed to be no spin 
order (^]. The theory of this paper is not literally ap- 
plicable to small-S" auantum systems, since no wavefunc- 
tion is possible in which every tetrahedron is simultane- 
ously a singlet. Nevertheless, it is claimed 48] that the 
diffraction from exact diagonalizations of the spin-1/2 
case agrees well with formulas such as (13.1411 . 

A pyrochlore Ising antiferromagnet, made into a quan- 
tum model by a small transverse ring exchange, has also 
defined the same coarse-grained field P(r), which in their 
theory is called an "electric field" ^60]. Out of their 
quantum- mechanical variables conjugate to P(r) they 
construct a "magnetic field"; the resulting U{1) gauge 
theory has gapless modes with |k| dispersion analogous 
to light waves, and correspondingly there are power law 
correlations (though with a different power law.) Ref. 45 
also studied a model with a polarization (they call P(r) 
the "magnetic" field) as a means to constructing quan- 
tum models with no long range order and fractionalized 
excitations. 

Recently, interesting phenomena have been observed 
in certain (electron) conductors containing a pyrochlore 
sublattice I69l|, which are speculated to be related to the 
frustration of this lattice. In particular, heavy fcrmion 
behavior is seen |69j | in the spinel LiV2 04, and unsual 
ferromagnetic behavior in pyrochlore Nd2Mo207 is as- 
cribed to a Berry phase acquired by the fermions in a 
spin background [71|. Perhaps the coarse-grained polar- 
ization field can help in modeling the long-wavelength 
behavior of these systems. 

There is also a speculation |70|] that the low- 
temperature state of ice itself (neglecting the dipole cou- 
plings beyond the nearest neighbor!) is dominated by 
proton tunneling with nontrivial Berry phases. Even 
though this model is (probably) not relevant to real ice, 
its exotic ground state is of interest in its own right and 
the P(r) field is likely to enter its decription. 

D. Diffraction singularities due to constraints 

Local constraints produce diffraction singularities in 
other systems, specifically in the so-called "transition 
state" of certain metal alloys [g^. In an FCC lattice, 
let (R) = ±1 represent two different chemical species. 
To model a state which has strong short-range order, as- 
sume a constraint resembling (|1.5|) . 

^s(R + u™)=0, (4.3) 

m 



where {R-(- u„j} are the 12 nearest-neighbor sites. After 
Fourier transforming, we obtain 

F(k)5(k) = 0, (4.4) 

where F{k) = exp(zk-Um). It follows from 14.411 this 
that the diffuse scattering (|s(k)p) is zero everywhere in 
reciprocal space, except that it diverges 62] along the 
two-dimensional surfaces defined by -F'(k) = 0. When 
these surfaces intersect the Ewald plane of an electron 
diffraction experiment, they produce well-known arcs ob- 
served in fee metal alloys near ordering transitions |62j] . 
Similar behavior is seen in the Na22Bai4CaN6 struc- 
ture, a triangular arrangement of rods each having an 
Ising degree of freedom with "antiferromagnetic" corre- 
lations J2\. Arcs arc also seen in quasicrystals, where 
they are ascribed to the constraints of tiling space > 

Eq. 1)4.4(1 is a sharper singularity than is found for the 
pyrochlore lattice in this paper. The fundamental reason 
is that the "spins" in (|4.3|l are on a primitive Bravais 
(FCC) lattice, so in a lattice of N cells there are N con- 
straints and the same number of variables. In contrast, 
in the pyrochlore problem analyzed in the present paper, 
there are four sites per primitive cell and only two con- 
straints. Consequently 1)4.4(1 gets replaced by the 4x2 
matrix equation ((3.51) . Just as 1(4.41) has its singular sur- 
faces at k values where the number of constraints, so ((3.5() 
has its pseudodipolar singularities at the points (recipro- 
cal lattice vectors) at which the two equations are linearly 
dependent and reduce to one equation. 

The above constraint-counting arguments have been 
phrased so as to make clear how they might be adapted 
to other systems. For example, the triangular Ising anti- 
ferromagnet (in zero field) is a highly frustrated system 
on a Bravais lattice, so one might naively expect stronger 
kinds of singularity in its diffuse scattering. However, 
in that case there is no equality constraint but instead 
Tliiea ~ whole approach breaks down. 
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FIG. 1: (a) Ising ground state on fragment of pyrochlore lat- 
tice; spins are shown on edges of a diamond lattice, (b). The 
corrsponding ice-model arrows between the diamond vertices. 
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2: Calculated structure factor in the plane normal to 
. (The figure is the same slice as shown in Refs. and 



